Temporal stimulated intersubband emission of photoexcited electrons 
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We have studied the transient evolution of electrons distributed over two levels in a wide quantum 
well, with the two levels below the optical phonon energy, after an ultrafast interband excitation and 
cascade emission of optical phonons. If electrons are distributed near the top of the passive region, 
a temporal negative absorption appears to be dominant in the intersubband response. This is due 
to the effective broadening of the upper level state under the optical phonon emission. We have 
then considered the amplification of the ground mode in a THz waveguide with a multiquantum well 
placed at the center of the cavity. A huge increase of the probe signal is obtained, which permits 
the temporal stimulated emission regime of the photoexcited electrons in the THz spectral region. 

PACS numbers: 73.63.Hs, 78.45. +h,78.47.+p 
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I. INTRODUCTION 

The temporal emission of THz radiation due to the 
coherent oscillations of electrons under ultrafast inter- 
band excitation have been studied during last decades 
(seei for review). The duration of the generated THz 
pulse is in the order of picoseconds due to effective relax- 
ation of the coherent response. Furthermore, the steady- 
state spontaneous and stimulated emissions from differ- 
ent semiconductor structures under an electric field pump 
have been demonstrated (see^ and Refs. therein). How- 
ever, the character of the temporal evolution of photoex- 
cited electrons during the non-coherent stage of the re- 
laxation is not completely understood yet. After the ini- 
tial photoexcitation and emission of the optical phonon 
cascade, which takes place in a picosecond time inter- 



val, which have been considered in details 



.1.3.4.5 



and 



Refs. therein), a non-equilibrium distribution appears 
in the passive region, with the energy Sex less than the 
optical phonon energy, hiOo^. A temporal evolution of 
this distribution takes place during a nanosecond time 
interval. This evolution is caused by the quasi-clastic 
scattering of electrons with acoustic phonons. Due to 
the partial inversion of this distribution, a set of pecu- 
liarities of the magnetotransport coefhcients, such as the 
total negative conductivity, or the negative cyclotron ab- 
sorption, appears^. 

For the wide quantum well (QW) case under consider- 
ation, with two levels in the passive region, the character 
of the evolution appears to be more complicate due to 
the interlevel scattering. The regime of negative absorp- 
tion is possible due to a more effective broadening of the 
absorption (dashed arrow in Fig. la) in comparison with 
the intersubband emission (solid arrow) . Such regime ap- 
pears because the absorption process involves the state 
in the active region for which the optical phonon emis- 
sion is allowed. Thus, a question arises about the tem- 
poral stimulated emission coming from a wide multiple 
quantum well (MQW) structure placed at the center of 
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FIG. 1: Scheme of the intersubband transitions for electrons 
with energies near the top of the passive region (the stim- 
ulated emission and absorption are shown by the solid and 
dashed arrows, respectively) (a), the geometry of the THz 
waveguide with the MQW structure in the center (6), and the 
initial distributions over the passive region for the cases A 
and B (c). 



THz waveguide, as it is shown in Fig. \h. In this pa- 
per we consider the amplification of the probe ground 
mode in the THz resonator caused by the temporal neg- 
ative absorption described above. We have found a huge 
amplification of the probe signal, so that it may be con- 
cluded that a temporal stimulated emission takes place 
in the THz waveguide with weak cavity losses and with 
the MQW structure placed at the center. 

Our study is based on the general kinetic equation for 
the distribution function in the conduction c-band, ]at^ 



written in the following form 
dfat 



dt 



Ga.t + Jif\at), 



(1) 



where Gat and J{f\at) are the photogeneration rate and 
collision integral for the c-band state a, respectively. The 
generation rate in (1) is given by*: 



Gat = wl ^ 5/(5a, (£a - £/) 



(2) 



where Wt is the form-factor of the excitation field with 
the pulse duration r^, the factor gi describes the relative 
contributions of different valence u-band states (heavy 
and light holes confined in the wide QW). The shape 
of the photoexcited distribution is given by the Gaus- 
sian function (5a(£^) = exp[-(£;/A)2]/(0FA) with the 
width A and centred at £■ = 0. The width of these 
function is determined by two different broadening pro- 
cesses: the spreading due to the energy-time uncertainty 
relation and the anisotropy of the valence band. These 
values are estimated by the energy values h/Tp and (iei, 
respectively, where /3 <C 1 is the anisotropy parameter. 
Under the optical phonon emission an additional broad- 
ening appears. Such a broadening is caused by the weak 
phonon dispersion, and characterized by the energy 6eopf 

Thus, the initial conditions in the passive region, af- 
ter the picosecond stage of the evolution due to optical 
phonon emission, are obtained from Eqs. (1,2) in the 
following manner: 






0)(£-£F')- 



(3) 



Here p^^ is the 2D density of states, the total broadening 
energy is determined as Aj^' ~ max(?i/rp, /3e/, 5eopt), 
and j = 1, 2 are the corresponding two levels in the pas- 
sive region. Below, we will concentrate solely on two 
models of the initial distribution: a narrow Gaussian 
peak (case A) and a flat distribution over the interval 
{HlOo — £21, f>^o) (case B), which are shifted from the 
boundary of the passive region, as Fig. Ic shows. Both 
distributions are relaxed between the sub-bands and to 
the bottom of the passive region. The temporal intersub- 
band response is determined by the interlevel redistribu- 
tion of the population. 

The analysis we will carry out next is divided in two 
sections. The temporal evolution of the distribution 
over the two subbands is described in Sec. II, includ- 
ing the consideration of the resonant intersubband re- 
sponse. The results for the transient amplification of a 
probe mode in the THz waveguide are given in Sec. Ill, 
with our conclusions presented in Sec. IV. The Appendix 
contains the microscopical evaluation of the broadening 
energy for the intersubband transitions. 



II. TEMPORAL EVOLUTION OF THE 
ELECTRONIC DISTRIBUTION 

We shall now turn to consider the temporal evolu- 
tion of the photoexcited electrons in the passive re- 
gion caused by the quasielastic scattering with acoustic 
phonons. Since the above-discussed initial distributions 
are isotropic over the 2D-planc, one has to consider the 
energy dependent distribution functions fi^2et governed 
by the system of kinetic equations: 



df: 



let 



at 

df2et 
dt 



Ji.f\^£t) - ^e-e-21 if let ~ he-e^it) , 
J{f\2et) + 1/,+,,, (/le+e.it - het) ■ (4) 



Moreover, these equations are written for the intervals 
(0, fiWo) and (£21 — hjjo, f)iOo) for the first and second sub- 
bands, respectively. The interlevel relaxation frequency 
is introduced here as 



27r 



p',q± 

x|(2|e-'^''-^|l) I' (2Ar« + l)<5 (£-£'), (5) 

where Cq is the matrix element for the bulk electron- 
phonon interaction with the acoustic phonons. The 
quasielastic collision integral in the jth subband is 
reached as 



Aim = I 



D 



df: 



jet 



3^ 



de 



V,ef: 



jeJjet 



(6) 



The diffusion and drift coefficients in (6), Dj^ and Vje, 
are determined as follows: 






7rP2 



2tt 



dtp 
2^ 



X \{j\e-''^'-\j)\ 



, 2n 

{2Nq + 1) ihw^)y2 

huJn 



(7) 



where Nq = [exp [hhoq/T) — 1] .A similar description 
of the quasi-elastic relaxation in a two-level system was 
discussed in Ref. I^- 

Since the elastic interlevel relaxation is the dominant 
process for the energy interval (e2i,fi^o)j one obtains 
/let ~ f2e-e2it — fet if t cxceeds the interlevel relax- 
ation time. Using the energy variable 1^ = e — hiOo we 
determine the distribution function f^+hujat = fit from 
the diffusion-drift equation for the two-level zone of the 
passive region, > ^ > £21 ~ hwo'- 



dt 



K^'^t-"^^^' 



(8) 



where D^ and V^ are the energy diffusion and the drift 
coefficients determined by Eq. (7). These coefficients 
are shown in Fig. 2 for a Ga^s-based QW which is 320 
A wide, with an interlevel distance £21 = 15 meV and 
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FIG. 2: Normalized diffusion and drift coefficients versus en- 
ergy. Solid and dotted curves: diffusion coefficient for T —4.2 
K and T =20 K respectively. Dashed line: drift coefficient 
which is independent of T. 
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FIG. 3: Evolution of the electron distribution for the case 
B. Solid, dashed, and dotted curves correspond to t —0 ns, 
1.6 ns, and 3.2 ns respectively. Vertical arrow indicates the 
position of £21. 



for two different temperatures T = 4.2 K and T = 20 
K. Thus, the normalization coefficients D^w^ = D^^q = 
3.91 X 10^ meVVs and 9.55 x 10^ meV^/s for T = 4.2 K 
and T = 20 K, respectively. Another coefficient, which is 
independent on the temperature, is Vuljo = V^=o = 5.22 
X 10^ meV/s. Note, that D/V ~ T for f = 20 K while a 
visible distinction of energy dependencies appears at T = 
4 K. We use below the energy-independent coefficients 
F = 3.5xl0^meV/s, andi:) = 2.2x10^ meVVs and 6.2x 
10^ meVVs for T = 4.2 K and T = 20 K, respectively. 
This approximation is valid within an accuracy better 
than 20% for the numerical parameters in the interval 
(hujo — S2i,fiuJo) considered below. Another important 



parameter is the interlevel relaxation time. 



which 



was assumed to be shorter than the time scales of Figs. 
3,4. Actually, for the above used parameters, i'~^ varies 
from 0.21 ns to 0.19 ns over the passive region. 

Equation (8) may be considered with the zero bound- 
ary conditions at ^ ^ ±00 if fkUo — £21 > £21- The 
initial distribution is located in the region ^ < 0, and 
electrons do not reach the bottom of the second level, 
i.e., £, > —£21- For the case A we use the initial con- 
dition: f^t=o = nexSAi^ - £.ex)l'^PmDi whcre riex is the 
total excited concentration, ^^^ ~ —9 meV is the exci- 
tation energy, and the half- width of the peak is equal to 
2.5 meV. Thus, the solution for Eq.(8) may be written 
as the moving peak 



fit 



2ft 



-SaA^ + VI-^,,) 



(9) 



fi2D 



with the time-dependent half-width A* — VA^ + ADt. 
For the case B, the broadening of the stepped distribu- 
tion is equal to 2.5 meV. The solution is given by the 
integral of the initial distribution multiplied by a similar 
to (2) form-factor. This solution is plotted in Fig. 3 for 
different times. 
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FIG. 4: Temporal evolution of the concentration for the dis- 
tributions A and B, at T =4.2 K and 20 K (solid and dashed 
curves, respectively) . 



In order to describe the temporal negative absorption, 
we need the concentration over the negative absorption 
region {hujo - 621, h^o) defined as rit = 2pg2D J^^^^ d^ht- 
The temporal evolution of nt, obtained with the solutions 
of Eq.(8) for the cases A and B, and the QW parameters 
detailed above, is shown in Fig. 4. 

The temporal evolution of the resonant absorption of 
the non-equilibrium electrons is described by the real part 
of the conductivity (see® and Appendix) : 



Re(TSe,t = 



e^Wi_?nt 



7 



£2l/?l 2 Vfe2+r2 fc2-K72 



(10) 



Here 5e — hw — 821 is the detuning energy, wb_l is the 
intersubband velocity matrix element, F and 7 are the 
broadening energy values due to the optical phonon emis- 
sion and the elastic scattering, which correspond to the 



absorption and stimulated emission processes. Using 
Eq.(All) and the parameters given above we make an 
estimate of F ^ 1.3(1.5) meV [for ^^^ — — 4(— 9) meV, re- 
spectively], so that r ^ 7 and one obtains Rease=o.t < 0. 
The contribution to Rease.t from electrons distributed 
over the region —£21 > 1^ > £21 ~ fii^^o is equal to zero 
due to the same broadening of the absorption and emis- 
sion processes. We have also supposed that electrons are 
absent from the region ^ < £21 — hcuo, where only the 
absorption due to transitions 1 ^ 2 is possible. 



III. TRANSIENT AMPLIFICATION 

We then consider the amplification of a probe THz 
mode in the ideal waveguide of width L due to the tem- 
poral negative absorption discussed above. The trans- 
verse electric field in the resonator E^^ exp(— itjt + ikx) 
is ruled by the wave equation: 






E' 



C2 dt 



= 



(11) 



with K^ = k"^ ~ euP' I (? ^ where e is the dielectric per- 
mittivity supposed to be uniform across the structure. 
The boundary conditions at 2; = ±L/2 takes the form 
-^zi±L/2 t ~ 0- ^* *^^ center of the resonator [z — 0), 
where the MQW is placed, one has to use: 



;iB_L 



^i* Lo = 0' 



dE: 



iBJ_ 



dz 



-iN- 



^Attuj 



-ase.tE'^iou (12) 



where N is the number of wells in the structure. The 
initial condition for the ground mode propagating along 
the resonator is: -Ef^^o — Ecos{ttz/L). 

Taking into account the slowness of the temporal evo- 
lution under the condition a = NTre'^\vs!,±\'^Lnex/jc^ ^ 1 
and restricting ourselves to the resonant case, Se — 0, 
we obtain the solution of Eqs. (11, 12) as E^^ = 
£t cos{t:z/L) where the time-dependent field is governed 
by: 



nt 



dt T*nex 



1 



4q:c"^ 

£0^21^^' 



(13) 



with the initial condition: £t=o = E. 

After averaging over the resonator, the Poynting vector 
along OX-direction is introduced as 



2ttu! 



dz\£,t? 



(14) 



and the temporal amplification is determined as the ratio 
of (14) to the Poynting vector at t = {): 



Vt 



St 
St=o 



exp(2r?^). (15) 



In Fig. 5 we plotted the logarithm of the amplification co- 
efficient versus time for a five-layers MQW with 7 = 0.2 







FIG. 5; Transient amplification rj^ for the same conditions as 
in Fig. 4. 



meV under the photoexcited concentration Uex — W^^ 
cm~^. We have also used e = 12.9 and L = 19 /im, cor- 
responding to the ground mode propagation along the 
THz waveguide. One can see the saturation of In 77^ when 
'T-t/'T-ea; approaches to zero (cf. Fig. 4). The maximal value 
of ln?7( is around 10^, i.e., a huge temporal amplification 
of the probe signal has been obtained. Thus, the stimu- 
lated regime of the emission should be realized in the THz 
waveguide with a weak damping (such waveguides have 
been studied recently, see Ref. |2l and references therein). 



IV. CONCLUSIONS 

In summary, we have considered a new non-coherent 
transient mechanism of the stimulated emission in the 
THz spectral region. The mechanism neither requires 
coherent response nor inverted distribution but it rather 
appears due to the different broadening of emission and 
absorption contributions in Eq. (10). 

The consideration performed here is based on several 
assumptions. Rather than using a microscopic calcu- 
lation of the photogeneration, we have used two mod- 
els for initial conditions {A and B) to obtain the tem- 
poral evolution of the distribution. We have also ne- 
glected the energy dependencies of 13^ and V^ in Eq. 
(8). In Eq. (10), which describes the intersubband 
response, we have neglected the Coulomb renormaliza- 
tion effect and taken into account that F 3> 7 as it is 
demonstrated in the Appendix. We also neglected non- 
parabolicity of the electron dispersion law and the ef- 
fect of the spatial confinement on phonons and electron- 
phonon interaction^". These approximations are gen- 
erally accepted for the GaAs/AlGaAs heterostructures 
with the parameters used here. In addition, instead of 
a stimulated emission effect due to the transient nega- 
tive absorption, we have calculated the amplification of 
a probe signal, and the role of the waveguide losses is not 
considered here. 



To conclude, we have found a huge amplification due 
to the temporal negative absorption under intersubband 
transitions of photoexcited electrons in a wide MQW 
structure placed at the center of a THz resonator. We ex- 
pect that the presented analysis motivates an experimen- 
tal treatment of the transient stimulated emission both 
for the system considered and for another heterostruc- 
tures with closely-spaced levels (stepped QW, tunnel- 
coupled structures, etc.). 
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APPENDIX A: BROADENING OF 
INTERSUBBAND TRANSITIONS 

Below we present the evaluation for the resonant in- 
tersubband conductivity given by Eq. (10). We consider 
the interaction of electrons with optical phonons at zero 
temperature and demonstrate that only the broadening 
of the 1^2 transition, which is responsible for the ab- 
sorption, takes place in the passive region while the 2— >1 
transition remains narrow (the broadening is determined 
by the weak elastic scattering). 

Within the framework of nonequilibrium diagram 
technique^'''^-:, one describes the transient response on the 
resonant probe field Es± exp(— iwi) based on the linear 
combination of Green's functions: 



Fi 



tlt2 



G 



tlt2 ' 



ARR 



de 



,ie(ti-t2) 



J- ti^ 



-5!Ftit2- 

Here T^t = 2f^t{G^"^ — Gf^) describes the temporal evo- 
lution of energy distribution f^t governed by the qua- 
siclassic kinetic equation (1), Q^^^'^ is the advanced 
{A) or retarded (i?) equilibrium Green's function, and 
6J^tit2 describes the linear response on the perturbation 
5h exp(— iwi). Using the e, i-representation we obtain the 

linearized equation for (5J^Eexp(— iwi) in the following 
form: 






(A2) 



where S^^''^ is the R- or A- self-energies and the integral 
contribution to this equation is given by 



Q 



-hujo ^ 



-iQr 



(A3) 



Here C, 



B(po) 
BQ 



is the matrix element of the Frohlich inter- 



action with the bulk phonon mode characterized by the 
wave vector Q = (q, q^^^ ) . 

For the case of the resonant intersubband excitation, 

\fujj — £21 1 ^ £21, only the component {6^^)21,1, = SF^p 
appears to be essential. Using the non-diagonal matrix 
element Sh2i — {ie/uj)EQ±vsi_ we transform Eqs. (A2) 
and (A3) into 



hiu- 621+ hp) SF^p - ^ Wpp'SFe^fiLOap 
p 

r ( riiiA r^RB. \ 

Je-hu/2 y-'e-hLj/2 ^ ^s-huj/2j^_^ 

f ( r^RA ri^iR 

Je+hui/2 \ye+huj/2 ^ '^e+hi^/2 



= 2—Eii±Vii 



22, p 



P 

(A4) 



where the broadening energy is given by 



\^e+n^/2) y^e+n^/2) 



ll,p. 



jQp' 

e — hhj /2 — hujo + «A 



|(2p|e»Q--|jp': 



e + hiij/2 — huo + iX 



(A5) 



and Wpp' is determined by a similar expression which is 
not essential below. The lower Eq. (A5) is written in the 
Born approximation and, from the formal point of view, 
A -^ +0. If we take into account the elastic scattering, A 
should be replaced by the elastic broadening energy. 

The Fourier component of current density due 
to the intersubband transitions under consideration 
is expressed through 6Fep according to j^^ ~ 
— iewi3j_ J dp J de6F^p/{2TThy. Here we have replaced the 
perturbation of the density matrix Spf using the relation: 
Sp^ — {~ih/2) linitj 2 5Ttit2- Performing the substitution 
5F^p — i2eEQ±VQ±ip /uj we obtain the conductivity : 



e^K'fiiP f dp 



{2TThy 



def. 



ep- 



(A6) 



The function ip is governed by the integral equation 

{hw - £21 + ilp) <y3ep - Yl ^PP"fie-huj^p 



f ( r^i,A r<^R 

Je-hLj/2 I ^e-hu/2 " '^e-hui/2 j , , 

\ / 11, p 

-Je+huj/2 \y'e+hu/2 ^ ^e+huj/2 



22, p 



(A7) 



which is obtained from Eq. (A4). 

Next we write the resonant conductivity (A6) for 
the case of the 5-like electron distribution /^ — 
(nf,x/2pQ2D)H£ex — £) with the excitation energy e^x lo- 
calized in the interval {hiOo — £ex, fiujo)- Due to the weak- 



ness of the integral term in Eq. ( A7) we obtain the con- 
ductivity in the form: 



a - ^>5£^n / de 



( ri^A 






huj - 621+ il' 



( f-i&A _ r-i&R 

\ye-hu/2 ^e-hu/2 



22, p 



hw-e2l +ijp \e^e,^~hu;/2 



P le— >£ea::+?it^/2 



(A8) 



The spectral densities in the numerators, 

i2X 

M,P (Sex - SjpY + A 



^R 



\pe-huj/2 ^e-huj/2 



, (A9) 



are written through the (5-functions if A ^ +0. As a re- 
sult, Eq.(A8) is transformed into Eq.(lO) with the elas- 
tic broadening energy in the second contribution while 
the optical phonon induced broadening appears in the 



first term. The corresponding broadening energy is de- 
termined as r = i?e7p|ep_^£^^,£^e^^+;ii^/2, or 



xS{£ex + £21 - ^^o - eip')|e _ 



(AlO) 



where V is the normalization volume and the effective di- 
electric permittivity, e* = {to — €oo)/£oioo, is introduced 
through the static and high-frequency dielectric permit- 
tivities, Co and Eoo- Performing the integration over the 
energy one obtains: 



e^hjJo 
4e* 



dqj_ J I 



d(p 






(All) 



9^ 
~^ex 



where the denominator contains the factor Qn, = {p^ 
Pg^ — 2pexpex cos ip) / fi^ with the characteristic momenta 
Pex = \/2meex and pex = ^2m[eex - (fi^o - £21)]- 



* Electronic address: ttvasko@yahoo. com| 
' Electronic address: ajhernan@ull.es 

^ J. Shah, Ultrafast Spectroscopy of Semiconductors and 

Semiconductor Nano structures (Springer, Berlin, 1999); F. 

Rossi and T. Kuhn, Rev. Mod. Phys. 74, 895 (2002). 
^ E. Gornick and R. Kerstling, Semicond. Semimet. 67 389 

(2001); V.N. Bondar, A.T. Dalakyan, L.E. Vorobev, D.A. 

Firsov, and V.N. Tulupenko, JETP Lett. 70, 265 (1999); 

M.A. Odnoblyudov, I.N. Yassievich, M.S. Kagan, and K.A. 

Chao, Phys. Rev. B 62, 15 291 (2000). 
^ H.Haug and A. -P. Jauho Quantum Kintetics in Transport 

and Optics of Semiconductors (Springer, Berlin, 1996). 

* F.T. Vasko and A. Kuznetsov, Electronic States and 
Optical Transitions in Semiconductor Heterostructures 
(Springer, New York, 1998). 

^ F. Prengel and E. SchoU, Phys. Rev. B 59, 5806 (1999); 

J. Schilp, T. Kuhn, and G. Mahler, Phys. Rev. B 50, 5435 

(1994). 
® S.E. Esipov and Y.B. Levinson, Adv. in Phys. 36 331 



(1987). 

■^ F.T. Vasko, JETP Lett., 79, 431 (2004). 

* A. Hernandez-Cabrera, P. Aceituno, and F.T. Vasko, Phys. 
Rev. B 67, 045304 (2003); S. Khan-ngern and LA. Larkin, 
Phys. Rev. B 64, 115313 (2001). 

^ B.S. Williams., S. Kumar, H. Gallebaut, Q. Hu, and J.L. 
Reno, Appl. Phys. Letters 83, 5142 (2003); H.N. Rutt, 
Z.J. Xin, and LA. Tan, J. Phys. D-Appl. Phys. 35,1907 
(2002); T.G. Ulmer, M. Hanna, B.R. Washburn, S.E. 
Ralph, A.J. Spring- Thorpe, IEEE J. Quantum Electron- 
ics 38, 19 (2002); M. Rochat, L. Ajih, H. Willenberg, J. 
Faist, H. Beere, G. Davies, E. Linfield, and D. Ritchie, 
Appl. Phys. Letters 81, 1381 (2002). 
^^ M.A. Stroscio and M. Dutta, Phonons in Nanostructures, 
(Cambridge University Press, Oxford, 2001); P. Kinsler, R. 
M. Kelsall, and P. Harrison, Physica B 263, 507 (1999). 
^^ E.M. Lifshitz and L.P. Pitaevskii, Physical Kinetics (Else- 
vier, Amsterdam, 1981). 



